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We relate short-range correlations in ultracold atomic Fermi gases to the entropy of the system 
over the entire temperature, T, vs. coupling strength, —l/kpa, plane. In the low temperature 
limit the entropy is dominated by phonon excitations and the correlations increase as r". In the 
BEC limit, we calculate a boson model within the Bogoliubov approximation to show explicitly 
how phonons enhance the fermion correlations. In the high temperature limit, we show from the 
virial expansion that the correlations decrease as 1/T. The correlations therefore reach a maximum 
at a finite temperature. We infer the general structure of the isentropes of the Fermi gas in the 
T,—l/kFa plane, and the temperature dependence of the correlations in the unitary, BEC, and 
BCS limits. Our results compare well with measurements of the correlations via photoassociation 
experiments at higher temperatures. 
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I. INTRODUCTION 

The thermodynamics of an ultracold atomic gas inter- 
acting via short range interactions is encoded in the short 
range two-body correlations of the particles. Such con- 
nections between short range correlations and physical 
properties of systems with short range interactions have 
been long recognized, e.g., Ref. and have been em- 
ployed recently to study the rf spectra of paired Fermi 
gases 0, H, I3- Tan has derived a number of remark- 
able relations between the properties of ultracold atomic 
gases, e.g., the ground state energy and their short range 
correlations 0] (see also [^). These links open new per- 
spectives to examine the many particle physics of systems 
with short range interactions. 

By using Feshbach resonances to tune the scattering 
length a, one can study atomic Fermi gases through 
the crossover from a weakly coupled BCS superfluid to 
a weakly coupled Bose-Einstein condensate (BEC) of 
molecules. In the strongly correlated unitarity regime, 
\a\ oo, the many-body correlations are highly non- 
trivial and one has to resort to Monte-Carlo calculations 
to obtain controlled results 0,1- At non-zero temper- 
atures T, information on the correlations is even more 
limited. We study here how knowledge of the free en- 
ergy and entropy can illuminate the nature of the finite 
temperature short range two particle correlations. The 
two-body correlation function in a two component Fermi 
gas has the general structure 



7,0,(r)0*(r'), (1) 



as one sees by regarding the correlation function as a 
Hermitean operator in r,r'. For r ^ d, the interparti- 
cle spacing, the functions essentially s-wavc Jastrow 
factors in the many-body wavefunction at short inter- 
particle distance, are determined by two-body physics, 
and outside the range of the potential, rg, at temper- 



atures of interest, are ~ sin(fc7' + 6)/r — sin 5 x(^)/''i 
where x{i") = 1 ^ ^/"^i (5 is the s-wave phase shift, and 
cot i5 = — 1/fca, where a is the s-wave scattering length. 
Thus 



r a 



(2) 



The correlation (or contact) strength C, which is deter- 
mined by many-body physics, can be measured directly 
in photoassociation experiments 0, [lol | . 

Here we examine the temperature dependence of the 
correlations by calculating C(T) from the finite tem- 
perature thermodynamics. Quite generally, the short 
range correlations are related to the free energy density 
/ = Tre~^/"^/f2 (where fl is the volume of the system) 
by i: 
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Differentiating with respect to the temperature, T, we 
find 
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m d df 



ds 
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thus linking the variation of C with temperature to the 
coupling constant dependence of the entropy density s. 
The connection is illustrated below in Fig. [ijwhich shows 
the temperature dependence of C, and the related Fig. 
[21 which shows the isentropes in the T, —l/kpa plane. 

We focus on the unitarity, BEC, and BCS limits and 
provide controlled results for C for both high and low T. 
Using the virial expansion, we show that the correlations 
decrease with C oc v? jT in the high temperature limit, 
where n is the density. From the relation ([4]) between the 
T dependence of C and the entropy density, we show that 
C(r) increases with temperature as for low T, owing 



2 



to the phonon contribution to the entropy. In the BEC 
hmit, we assume a bosonic model to relate the bosonic 
and fermionic correlations; a calculation within the Bo- 
goliubov approximation illustrates how phonons enhance 
the fermion correlations at low temperature. Combining 
the behavior of C in the high and low T limits, we in- 
fer that C reaches a maximum value at a temperature 
Tmax > 0. This intriguing behavior of C{T) is detectable 
in photoassociation experiments near unitarity. 



II. CORRELATION FUNCTION AND FREE 
ENERGY 



We consider a homogenous two-component Fermi gas 
described by the Hamiltonian [h = 1 throughout) 



+ Jd'rJ d'r'V{\r-r'\)i:l{r)4{r')Mr')Mr), 

(5) 



where the ipi are the field operators for fermions in in- 
ternal state We take the range of the interaction 
V(r) to be short compared to the interparticle spacing d. 

We first present an elementary derivation of the rela- 
tion ^ between the two-body correlation function and 
the free energy density, using the device of scaling the 
potential V{r) by a factor A in dS]); V{r) XV{r), let- 
ting A —> 1 in the end. The correlation function for given 
A has the structure of ^ with x, C and a functions of A, 
For r <^ d, the function x\ satisfies the low energy limit 
of the two-body Schrodingcr equation [1, [ll[ 



strength is 



1^ 

m dr^ 



+ AF(r) xaW^O. 
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At short distances, with r > tq, Xa('') = ^ ^ f/o-Xi for 
convenience, we choose the normalization constant to be 
unity. Differentiating ^ with respect to A, multiplying 
by xa, integrating with respect to r from to rc, and 
eliminating the term J drx\Vdx\/dX using we find 



df 9/ dX _ 1_ /dH_\ _dX_ 
da^ ~dXda^ ^ n \ 9A / da^^ 
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and setting A = 1, we arrive at 



Combining ^ and 
the desired relation 

We review the application of ([3]) to a homogeneous 
equally populated two-component Fermi gas of volume 
ri, at T = The ground state energy per parti- 

cle E/2N can be written as (3/5) (1 + P)Ef, where (3 is 
a function of the dimensionless quantity ^ = —1/kpa. 
The Fermi momentum kp is defined in terms of the sin- 
gle component density n = N/Q by n = kp/dTv"^, and 
the Fermi energy by Ep = kp/2m. Equation ([3]) im- 
plies that C{T = 0) = {kj,/40TT^)dl3/d^. The function (3 
has been calculated by quantum Monte Carlo simulations 
near unitarity (1/a — *■ 0); at unitarity dp/dS, « 0.9 0, a 
result which agrees well with the value C = 2.7{kp/36Tr^) 
extracted directly from the correlation function calcu- 
lated by quantum Monte Carlo methods In the BCS 
limit, —1/kpa 3> 1, a perturbative expansion in kpa 
gives [H, M H 



— — X^/F 
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10, 
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(9) 



Corrections to E due to pairing are exponentially small 
in this limit. Keeping the leading term yields C = a^n? . 
In the BEC limit, 1/kpa S> 1, fermions form molecules 
and [13, m 



-.{kpar, 



\3/2 



(10) 



with Ei, = — 1/ma^ the molecular binding energy, and 
am the scattering length between molecules. The second 
term is the Hartree-Fock mean field energ y an d the third 
the nonperturbative Lee-Yarig correction [l3!|. Few body 
calculations give = 0.6a [1^. In the BEC limit, C = 
n/2Tia to leading order. The divergent behavior, as a ^ 
0, arises from the normalization of the molecular wave 
function, x/^ ~ e~'"/''/ra^/^. 



dX 



1 



drV{r)xi{r)—^= 



da^ 



m 



(7) 



where we use vanishing of the s-wave function xa at the 
origin, and extend the upper bound of the integral from 
Tc to infinity since V{r) is nonzero only for r < tq. On 
the other hand, from the Feynman-HcUmann theorem, 
the variation of the free energy density with coupling 



III. LOW TEMPERATURE 

We now consider the temperature dependence of C for 
T well below the supcrfluid transition temperature Tc, in 
which regime phonon excitations dominate. The contri- 
bution of the phonons to the entropy density is 



5pho 
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where the zero temperature sound velocity Cs is given by 
mc^ = ndjji/dn, with fi is the chemical potential. The 
change of C to leading order for T <^ Tp = Ep (fcs = 1) 
is thus 



SC = C(T) - C(0) 



irm 
120 



T 



(12) 



The zero temperature sound velocity appears to in- 
crease monotonically from the BEC side to the BCS 
side 0, [3, [i3| , so that the short range pair correlations 
as parametrized via C increase from zero temperature as 
At T = 0, At = (1 + /3 - so that 



1 



/V3 



= 1 + /3 - -/3'^ + T7r/3"^^ (13) 
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with f3' = df3/d£,, etc. Recent measurements of Cg are in 
good agreement with combined with Monte Carlo 
results for /3 around unitarity [2]| . 

Equation ([9]) implies that in the BCS limit, = 
(z;^/3)(l + 2kFa/TT), and thus at low T, 



c(r) = 



160 



TV 



(14) 



For Tc < T Tp, the weakly attractive Fermi gas 
can be described by Landau Fermi liquid theory. The 
entropy density is given by [2^ s/ = ^m*kpT with 
m* = TO[l + 8(71og2- l)(A:i.a)Vl57r2] the effective mass 
of the quasiparticles [2y|. From we then have the 
temperature dependence 



2„2 



SC{T)/n^a 



-^{7log2- l)kFaf^) >0. (15) 



In the B EC limit, the Bogoliubov sound speed is Cg = 
\/ ngm/J^i with gm = '^'Ka„jM and M ~ 2m (23|. Thus 



1 



da 

(4867ri3)i/2 
40(fcFa)5/2 
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At unitarity, 



80(l + /3)5/2 



T 



IV. CORRELATIONS IN THE BEC LIMIT 

The universal behavior 6C ^ T"* as T — > shows that 
thermally excited phonons enhance the fermion pair cor- 
relations. In the BEC limit, the enhancement of the cor- 
relations can be understood directly in terms of a gas of 
molecules interacting with a short range potential with a 



scattering length > 0. The boson correlation function 
is similar in structure to that for fermions: 



lim(0t(r)0t(o)^(O)c/.(r)) = C„, (-- — 



where the 
fermions. 



(18) 



(r) are the bosonic field operators. As for 



Cn 



M df„ 

271 dttm 



(19) 



with /,„ the free energy density for the bosonic molecules. 
The factor of 2 difference from ([3]) is due to the factor 
1/2 in the interaction energy for identical bosons written 
in terms of the (p. 

Since the low temperature physics described in terms of 
the fundamental fermions or bosonic molecules must be 
the same, the leading temperature variation 5f ~ f{T)~ 
/(O) is equal to Sf,n = fm{T) - /„(0). Therefore 



5C 



m d{Sf) 
4^ 9a-i 



m d{5fm) 9a„ 



47r da 
m da~^ 
2M da- 



da~ 

SCm, 



(20) 



with SCm ~ Cm{T) — Cm(0). We calculate the boson 
correlation function within the Bogoliubov approxima- 
tion using a pseudopotential Vm{r) = gmS[Y)- We divide 
the field operator 0(r) = 00 + 54>{v) into a condensate 
part 4>o and a fluctuation Scj), assumed to be small. To 
second order in the fluctuations, at small r, 

(0t(r)0t(o)0(O)0(r)) 

= nl + 2no[2(<S0t (0)50(0)) + (<S0t(i.)50t(o))], (21) 



where 



lim((50^(r)(50^(O)) 



(22) 
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^ yukVk p 1- 2ukVk{alak) 



(2 



(23) 



(17) here vl = {^k/Ek - l)/2, ^ = {^k/Ek + l)/2, ^fe - 
k"^ /2M -t- gmriQ, and the ak annihilate phonons with 



energy dispersion Ek = V^fc ~ {gmno)'^ ■ The 1/r 
divergence in {5(j)^ {r)5(j)^ {Q)) corresponds to the cross 
term in (fT8|) . The ~ l/r^ in (fTS]) appears from the 
higher order term limr^o(<5(/)1'(r)(50t(O)) ((50(r),50(O)) = 
lim^^o |(<5'/'(r)<5'/>(0)) P . . ., not included in ([^1]) . 

The lowest order Bogoliubov approximation produces 
the correlation function to leading order in (na^)^/^. 
In (HH), the r independent term is reliable to order 
(rifl^J^/^, the term proportional to l/'n}/'^r to order 
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(na^j)^/'^, and the term proportional to l/(n^/'^r)^ to or- 
der (na^)^/^. To derive the fuU Cm(l - Om/r)^ struc- 
ture to higher order in (na^)^/^ requires calculating 
terms beyond the simple ones in (|2ip . a task we de- 
fer. To extract C,„ from the r independent term in 
(f2T|) we write no = n - {S(j){OY S<l){0)) . At T ^ 0, 

1 /2 

Cm/ {namY = 1 + (64/3) (naf^j/7r) , which agrees with 
the calculation of Cm from p9|) using the ground state 
energy of weakly interacting bosons [cf . (fTO|) ] . The ther- 
mally induced change (5 Cm is 



2nal 



(Pk 



fc2 



(27r)3 2^;^ e^'-Z'T - 1 



60 c° 



where we consider only the phonon contribution {E^. ~ 
Cgk). Substituting into (PO)) . we obtain the same 
result as in p6)) . 

The dependence can be understood in terms of the 
effect of thermal sound waves on the eigenstates, 0i(r), of 
the two particle density matrix ([T]). At distances beyond 
the interparticle spacing, Eq. ^ for the </> (with A = 1) 
contains potential terms from the mean field gmf^ and its 
fluctuations. At fixed normalization, the latter change 
the magnitude of 4'i{r) at short distances by terms 
{{SnY) ~ T^, which translates into the dependence 
of C at low temperature. 



V. HIGH TEMPERATURE 

In the high temperature limit, T ^ Tp, where the fu- 
gacity z = e^/^ is small, the free energy can be calculated 
via a virial expansion. To second order in z, the parti- 
tion function for a two-component Fermi gas is given by 

[Hill 



where 



log(Z/Zo) = -^z%2 



\El\/T ^ f°° dk d(5(fc) ^„fc2/^^ 

TT dk 



(25) 



(26) 



is the second virial coefficient. Here the El are the bound 
energies of the two-body attractive interaction V{r), 5{k) 
is the s-wave scattering phase shift, and A = (27r/mT)^/^ 
is the thermal wavelength. In the BEC regime, we con- 
sider only the relevant bound state with Eb = — 1/ma^. 
In the regime T <^ l/mrp, the values of k entering 
(p6)) are ^ I/t'Oi so that cot(5(fc) = —1/ka; hence 
dS{k)/dk = -a/ ((fca)2 + 1), and 




'dt 



(27) 



a dimensionless function of A/a. Thus 

mkT d log Z /tt 2 x 2 ^^^2 



C 



Ann da- 



d{\/ay 



(28) 



where we replace z by the number density n = 
{T/n){d\ogZ/dii) = (z/A3) (1 + 23/2z&2). 

For T ^ Ta = 1/ 'ma^ (a 3> A) , the system is a weakly 
interacting Fermi gas for all coupling strengths. From 
([271) . db2/d{X/a) = l/v/27r, and thus 



c(r) 



2 \2 



77,2 A 



2n2 
mT 



(29) 



Equation (|29|) combined with p2)) leads to the important 
result that the correlation C has a maximum at a nonzero 
temperature Tmax- 

In the BCS limit, for Tf<^T <$: Ta, 62 = -a/V2A and 
C(T) = Ti^a^; this is simply the Hartree term in the low T 
result (|14p persisting to higher temperatures. From p4p . 
(fT5)) and we expect that Tmaa; '^-^ TV and that the 

relative size of the maximum, [C{T,nax) ~ C'(0)]/C(0) ~ 
|fci?a|, is small. 

In the BEC limit, when the temperature is lowered to 
T ^ Ta {X ^ a), the bound state starts to dominate in 
([27)) and &2 ^ cxp{X'^ /2TTa^). The correlation C increases 
from ~ n^A^ to n^A"^ exp(A2/27ra2)/a. The virial ex- 
pansion for the Fermi gas breaks down when 62 z ^ 1, 
which happens at T ~ Ta/ \og{l/kpa). At this point 
the number of the bosonic molecules becomes significant. 
Within the regime Tp <^ T <^ Ta/logil/kpa), a virial 
expansion can be carried out for the bosonic molecules. 
The calculation proceeds as for the fermions except that 
there is no bound state for the molecules and one has 
to add the binding energy —n/ma^ to the free energy 
density, yielding 



c(r) 



n 
2'na 



(kpa)^ f da 



da 



(30) 



The second term in ([50]) comes from the mean field en- 
ergy for a normal weakly interacting Bose gas, which is 
twice that of a condensed Bose gas [cf. ([TS]) ]. From (|16|) . 

and ([501) . "^"^ expect that C reaches a maximum at 
Tmax - ~ Tf, and that [C{Tmax) - C(0)]/C(0) < 1. 

In the unitarity limit, ([T7)) and ([^^ indicate that C has 
a maximum at Tmaa ~ Tp with [C(Tma2;)-C(0)]/C(0) ~ 
I. The pronounced maximum in the correlations in the 
strong coupling limit indicates the importance of phonons 
there. 

Figuredlsketches C{T) in the BCS, unitarity, and BEC 
regimes. The solid curves are interpolations between the 
low T results ([n| - ([T7)) and the virial expansion results 
((28l)-((3Ql) for T :$> Tp. Note that C is continuous at 
since the superfluid transition is second order. 

The link between the entropy and the short range cor- 
relations through ([U is illustrated in Fig. where 
we sketch the isentropes in the BEC-BCS crossover 
regime, as well as the superfluid transition temperature 
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FIG. 1: (color online) The correlation strength C{T) in the 
BCS (a), unitary (b), and BEG (c) limits. The dashed lines 
at high T are the virial results, and the dashed line at low T 
in (b) is from The maximum in C{T) for T = T,nax in 

the BGS and BEG limits is too small to be visible here. 



T/Tf 
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FIG. 2: (color online) Isentropes of a balanced two-component 
Fermi gas, from which one can infer the temperature depen- 
dence of C(T). 



Tc [T1,[T1,[23|. It follows from (g]) and the virial expan- 
sion result ([M)) that the isentropes has a negative slope 
in the T, —1/kpa plane for high T. Likewise, the increase 
in C(r) given by P^ - P?)) means that the isentropes has 
a postive slope for low T. The correlation strength C(T) 
has an extremum where the slope of the isentropes as a 
function of —1/kpa vanishes. Additional information is 
obtained from a recent experiment, where a Fermi gas 
was prepared in the BCS limit with initial energy Ei and 
entropy Si and then adiabatically tuned to the unitary 
point where the final energy Ef was measured |2l[. The 
temperature was deduced from T = dE/dS. Within the 
experimental regime the final temperatures Tf at unitar- 
ity were generally higher than the initial temperatures 
Ti in the BCS limit [21|. However, for the isentrope end- 
ing at Tf w 0.2Tp, Ti w Tf, which indicates that this 
isentrope (as well as neighboring ones) first bends up- 
wards and then downwards from the unitarity region to 
the BCS limit. Well below the BCS transition, where 
phonons dominate the entropy, the isentropes have posi- 
tive slope, not visible on the scale of Fig. ^ . 



VI. PHOTOASSOCIATION EXPERIMENT 

Photoassociation experiments provide a direct measure 
of the correlation strength [9l.[l0j. In a recent experiment, 
^Li atoms were trapped in the lowest two hyperfine states 
|1) and 1 2). Then the bare closed channel spin-singlet 
molecular state \tpv=3s) associated with the 834G Fesh- 
bach resonance in the 1-2 channel was excited by a laser 
field Ei to a spin-singlet molecular state lipv'^es), with 
linewidth 7 = (27r)11.7 MHz. The excited molecules were 
lost from the trap and the remaining atoms counted for 



various durations of the laser pulse i9j. By Fermi's golden 
rule, the local loss rate of the number of atoms, 2n, in 
the trap is 



r 



dt 
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(31) 



where = (■(/' „' =68 1 El • d|-!/)„=38) is the Rabi frequency, 
d is the atomic dipole operator, and Ub the local density 
of molecules in the closed channel. The density of closed 
channel molecules is related to the correlation function 
between the states |1) and |2) in the open channel by [l3| 



rib 



1 



47raj,gC 
mjjLb^B \abg 



1 



(32) 



where abg is the background scattering length, the molec- 
ular magnetic moment /if, is twice the Bohr magneton /i^, 
and AS is the width of the Feshbach resonance. 

In a homogeneous gas in the high temperature limit, 
we find from dST]), and ([32]) that F - n'^/T. In 

a trap, the local density approximation gives n(r) = 
iV((I'/Ar)^exp[-/3^^(r)] with V{v) = {m/2){ujlx^ + 
t^yy"^ + uj\z^)l2 and cD = i^uOx^yU^zY^"^ . The average 
of C over the trap, and thus the average rate F, is 
cx <D' 



''N'^/T^/'^, and Eq. jH]) has the solution 



m 

N{Q) 



1 I ^RNb. 



(33) 



Analogous expressions for the T — trap-averaged loss 
rate were given in p^ . 

In Fig. [HI we compare the high T result ([33]) with the 
experimental data obtained on the BCS side of the res- 
onance Q. Initially the temperature is below Tf, where 
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((3T|) underestimates the loss; however, owing to depletion 
the Fermi temperature falls below T, and ([55]) provides 
reasonable agreement with the data - with no fitting pa- 
rameters. 
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FIG. 3; (color online) Loss of "^Li atoms in the trap through 
photoassociation vs. the probe duration at S = 865 G. The 
solid curve is the numerical evaluation of (|33|l . The circles are 
the T = 0.75Tf experimental data from Ref. [^. 
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